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ABSTRACT 


The  dissociation  of  carbon  monoxide  diluted  with  noble  gases  has 
been  studied  in  a  shock  tube.  The  reaction  is  complex  and  appears  to  be 
a  chain  with  as  the  important  intermediate.  At  temperatures  around 
8000°K,  the  reaction  rates  found  to  be  consistent  with  the  data  are 
O  +  C+  M-CO  +  M,  kr  =  2x  iO-34  cm6  sec-1;  C  +  C  +  M  -* 

C2  +  M,  kr  =  3  x  10"33  cm6  sec"1;  C2  +  O  -*CO  +  C  ,  k  =  6  x  10"10 
cm3  sec  1 ;  and  02  +  C  -*  CO  +  0,  k  =  5  x  10  111  cm3  sec  1 . 
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INTRODUCTION 


The  high  temperature  chemistry  of  carbon  is  of  interest  in  several 
fields  such  as  combustion  and  high  speed  entry  into  some  planetary  atmo¬ 
spheres,  but  little  quantitative  work  has  been  reported.  Since  CO  is  the 
simplest  starting  material  which  is  easily  handled,  it  was  chosen  in  this 
work  to  determine  its  suitability  as  a  controlled  source  of  carbon  atoms, 
and  this  paper  will  report  a  study  of  its  dissociation. 

With  a  binding  energy  of  11.  1  ev,  it  is  necessary  to  provide  tem¬ 
peratures  greater  than  5000°K  to  cause  appreciable  dissociation  and  the 
shock  tube  technique  was  used  as  the  experimental  method.  The  dissociation 
reactions  of  most  of  the  permanent  diatomic  molecules  have  now  been  studied 
in  shock  tubes.  Under  favorable  conditions  rate  constants  estimated  to  be 

within  20%  *  of  the  correct  value  have  been  obtained.  It  is  rate  that  reported 

2 

rates  differ  by  as  much  as  an  order  of  magnitude. 

For  systems  consisting  of  a  highly  dilute  solution  of  the  reactant  in 
inert  gas  theoretical  calculations  are  now  available  and  the  results  tend  to 

3 

be  in  good  (a  factor  of  2)  agreement  with  experiments.  It  is,  however, 
necessary  to  know  the  reaction  mechanism  and  while  this  is  a  straightforward 
matter  for  homonuclear  molecules,  heteronuclear  molecules  may  have  much 
more  complicated  reaction  schemes.  Such  complications  have  led  in  the  past 
to  uncertain  analyses  in  systems  such  as  NO  and  have  necessitated  in  more 
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recent  experiments  very  careful  data  acquisition  and  interpretation. 


EXPERIMENTAL 


A  1.  5  inch  diameter  shock  tube  was  used  and  measurements  were 
made  transversely  behind  the  incident  shock  (about  11  feet  from  the 
diaphragm).  A  dump  tank  several  feet  from  the  observation  station  pre¬ 
vented  wave  reflections.  The  tube  was  made  of  pyrex  g’ass  and  radiant 
emission  measurements  were  made  in  the  infrared,  visible  and  ultra¬ 
violet. 

The  infrared  measurements  were  made  through  a  calcium  fluoride 
window.  This  was  contoured  to  fit  flush  with  the  inside  of  the  tube  and 
cemented  in  with  an  epoxy  resin.  A  16  inch  focal  length  gold  coated  con¬ 
cave  mirror  was  used  slightly  off  axis  to  collect  and  focus  the  light  onto  a 
gold  doped  germanium  detector.  The  CO  fundamental  at  4,  6  fi  was  isolated 
with  an  interference  filter  (Fig.  1).  The  image  magnification  was  unity 
and  the  entrance  aperture  to  the  detector  2  mm  square.  With  f/l 5  or 
smaller  optical  apertures  the  spatial  resolution  was  about  4  mm.  The 
electronic  rise  time  (10-90%)  was  about  1  fj.s. 

Time  integrated  spectra  were  taken  in  the  visible  and  ultraviolet, 
with  a  small  quartz  spectrograph.  The  shutter  was  opened  before  the 
arrival  of  the  shock  wave  and  closed  later.  The  spectra  were  taken  through 
a  quartz  window  or  the  calcium  fluoride  window.  It  was  generally  necessary 
to  superpose  several  shocks  on  one  plate.  The  spectra  ohowed  strong 
Swan  bands,  CN  violet  and  atomic  carbon  radiation  and  as  a  result  these 
wavelengths  were  also  monitored. 
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Fig.  1  CO  emissivity  at  T  =  5000  K  and  IR  filter  transmission 
characteristics. 
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FILTER  TRANSMISSION 


In  the  visible,  the  was  followed  by  a  monochromator  plus 
photo  multiplier  set  to  receive  5200-5088$.  The  CN  was  followed  through 
an  interference  filter  at  4216$.  In  both  cases  the  light  was  viewed  directly 
through  the  glass  walls  of  the  shock  tube. 

The  atomic  carbon  line  at  2478.$  was  monitored  with  an  RCA  6903 
tube  mounted  in  the  focal  plane  of  the  small  quartz  spectrograph.  The 
entrance  slit  was  pushed  as  close  as  possible  to  the  window  and  the  optical 
axis  was  perpendicular  to  the  tube.  The  entrance  optics  were  masked  down 
to  f/30  giving  a  spatial  resolution  of  about  2  mm. 

Velocity  me"  surements,  estimated  good  to  l/2 %  were  made  from 
the  response  of  thin  film  heat  gauges.  These  were  painted  onto  plugs  cut 
from  the  wall  of  a  spare  piece  of  tube.  The  plugs  were  placed  flush  with  the 
inside  wall  and  sealed  with  epoxy  cement.  The  differentiated  signals  were 
displayed  on  a  raster  trace  on  an  osiclloscope. 

Most  of  the  experiments  were  made  by  bursting  the  diaphragm  with 
cold  hydrogen  at  pressures  up  to  1400  psi.  Some  experiments  were  made 
with  a  hydrogen-oxygen  driver  diluted  with  helium  and  ignited  by  an  exploding 
wire  strung  down  the  length  of  the  driver. 

Gas  mixtures  were  made  from  Matheson  high  purity  argon  or  neon 
and  CP  grade  carbon  monoxide  by  partial  pressure,  in  a  mixing  tank,  several 
hours  elaspsing  before  any  mixture  was  used.  An  early  tank  of  CO  had  an 
impurity,  which  was  never  identified,  with  a  strong  infrared  spectrum  and 
which  gave  anomalous  spikes  of  C^  at  th  shock  front.  The  CO  used  in  the 
experiments  reported  here  did  not  show  this  behavior.  The  only  infrared 
spectra  were  the  CO  fundamental  and  overtone.  The  main  impurity  was 
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nitrogen  (460  ppm):  Carbon  dioxide  (350  ppm)  and  oxygen  (150  ppm)  were 
the  only  other  significant  impurities. 

Initial  pressures  were  generally  a  few  torr  and  shock  speeds  up  to 
5mm/jisec  were  used  giving  temperatures  from  about  5000-10,  000°K. 

RESULTS 

The  spectrograms  obtained  with  the  small  quartz  instrument  showed 
the  visible  region  to  be  dominated  by  Swan  and  CN  violet  bands.  In  the 
ultraviolet  the  Mulliken  bands,  the  atomic  carbon  line  at  2478.$  and  two 
band  systems  of  CO  were  identified,  the  fourth  positive  bands  in  the  region 
2000-2800$  and  the  (0,  1)  band  of  the  Cameron  system.  Other  systems  of 
CO  such  as  the  triplet  and  3A  and  Asundi  were  sought  but  never  positively 
identified. 

Calibration  of  Infrared  System 

The  infrared  system  was  sensitive  enough  to  be  used  with  mixtures 
diluted  to  0.  5%  CO  +  99.  5%  Ar  or  Ne  and  initial  pressures  of  5  torr.  At 
lower  CO  concentrations  the  signal  to  noise  ratio  was  too  low  for  meaningful 
data  to  be  obtained. 

Calculations  showed  that  even  at  the  greatest  CO  density  used  the 
intensity  of  infrared  radiation  at  temperatures  greater  than  4000°K  was 
proportional  to  the  temperature,  and  the  number  of  molecules,  i.  e.  ,  the 
line  emissivity,  was  «  unity.  The  detector  and  recording  system  should 
have  a  linear  response.  An  overall  check  was  made,  and  the  results  are 
shown  in  Fig.  2  where  the  results  encompass  a  temperature  range  from 
4000-6000°K.  (In  the  analysis  of  runs  in  which  reactions  occurred,  only 
relative  measurements  were  made.  )  The  experiments  were  made  well  in 
the  linear  part  of  the  curve  (p^AT.  10"  ^  <  5  where  p^  is  initialpressure 
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P|<£A  T  (  TORR  x  FRACTION  OF  COx  Ax  T  xIO"4) 

Fig.  2  IR  detector  system  output.  The  abscissa  is  in  units  of  torr 
x  fraction  of  CO  x  density  jump  x  T  x  10“  . 
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in  torr,  0  is  fraction  of  CO  x  100,  A  is  the  density  jump,  T  is  the  tem¬ 
perature.  )  Figure  2  correlates  data  up  to  10%  CO  in  Ar.  It  was  found  that 
pure  CO  runs  did  not  scale  in  this  way,  being  apparently  self-absorbed. 

The  addition  of  argon  pressure  broadened  the  lines  enough  to  make  them 
optically  thin. 

Vibrational  Relaxation 

The  data  of  Hooker  and  Millikan  ^  show  that  at  temperatures  above 
about  5000°K  in  the  present  experiments  the  vibrational  relaxation  time 
will  be  less  than  the  rise  time  of  the  experiment  and  will  be  unresolved. 

For  a  given  run  and  under  conditions  of  constant  pressure  behind  the  shock 
and  large  dilution,  the  infrared  signal  will  be  proportional  to  the  fraction 
of  original  CO  present. 

At  low  temperatures  it  was  found  possible  to  follow  the  CO  vibrational 
relaxation  and  the  results  were  in  good  agreement  with  the  results  of  Hooker 
and  Millikan.  At  higher  temperatures  the  rise  time  became  limited  by  the 
apparatus  and  at  still  higher  temperatures  (around  6000°K)  the  CO  concen¬ 
tration  began  to  drop  off  behind  the  shock  front. 

Incubation  Times 

The  dissociation  of  the  CO  as  observed  did  not  occur  directly  be¬ 
hind  the  shock  front  as  expected  but  after  a  noticeable  incubation  time. 

With  increasing  temperature  and  pressure  the  incubation  time  became 
shorter.  The  incubation  times  are  shown  in  Fig.  3.  Most  of  the  data  on 
these  incubation  times  were  obtained  with  argon  mixtures  and  frequently  no 
measurable  slope  was  displayed  during  the  incubation.  Some  runs  made 
with  neon  did  show  an  initial  decrease  in  the  infrared  trace  during  the 
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Fig.  3a  1%  CO  in  Ar,  =  10  torr,  Us  =  2.  47  mm/ps.  Lower  trace, 

IR  5  p  ;  Upper  trace,  visible  0.  5  p  :  time  base  20  p  s/div. 

b  2%  CO  in  Ar  pj  =  3.  1  torr,  Us  =  2.  97  mm/ps.  Lower  trace, 

IR  5p  ;  Upper  trace  visible  0.  5p,  time  base  lOps/div. 

c  10%  CO  in  Ar  p^  =  20  torr,  U3  =  3.  07  mm/ps.  Lower  trace, 
IR  5p  ;  Upper  trace,  visible,  0.  5 fj,  ;  time  base  10  p s/div. 

d  The  same  run  as  (b). 

Trace,  UV  2478^  ,  time  base  10  ps/div. 

e  1%  CO  in  Ar  p^  =  10  torr,  Us  =  3.  02  mm/ps.  Lower  trace, 
visible  0.  5 p  ;  Upper  trace  IR  5 p,  time  base  20ps/div. 


incubation  time.  With  the  incubation  times  as  defined  in  Fig.  3,  Fig.  4 
was  constructed  to  correlate  the  results. 

Concurrent  observations  on  the  atomic  carbon  line  showed  a  slow 
rise  followed  by  a  faster  one,  Fig.  3d.  The  incubation  times  data  obtained 
from  such  traces  are  given  in  Fig.  4c.  The  agreement  obtained  in  the  same 
experiment  of  incubation  times  from  UV  and  IR  measurements  was  never 
perfect  but  always  within  a  factor  of  two,  which  is  the  experimental  scatter. 

The  2478$  ob  servations  were  more  sensitive  at  low  CO  concentrations  and 
enabled  incubation  data  to  be  obtained  at  dilutions  of  0.  02%  CO  in  Ar, 

The  excitation  of  the  2478$  line  of  carbon  in  shock  tubes  is  not  fully 
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understood  and  the  subject  of  another  report.  While  the  emission  level 

depends  on  a  number  of  factors  it  is  expected  to  scale  in  these  experiments 

with  concentration  and  temperature. 

C  Bands 

_2 _ 

The  C^  Swan  bands  were  observed  to  rise,  pass  through  a  maximum 
which  always  occurred  later  than  the  end  of  the  incubation  time,  and  then 
fall  closely  following  the  CO  trace.  There  was  often  a  small  toe  at  the 
shock  front  in  the  C^  history  indicating  a  rate  of  production  proportional 
to  time.  In  quite  a  number  of  runs,  however,  the  toe  was  not  very  discernable. 

The  CN  violet  radiation  history  was  quite  similar  to  that  of  the  C^ 

Swan  bands.  It  was  only  studied  as  an  adjunct  to  some  experiments  where 
nitrogen  was  deliberately  added  as  an  impurity. 

Reaction  Mechanism 

The  presence  of  C^  shows  that  the  reaction  is  complex  and  the 
shuffle  reactions  may  be  important.  Unfortunately  oxygen  does  not  have 
a  suitable  spectral  region  for  study  which  will  not  be  dominated  by  other 
radiation  under  the  conditions  of  these  experiments,  and  so  it  was  not  studied. 
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Fig.  4a  Incubation  times  from  infrared  experiments  plotted  as  function 
of  rp,  (0)1/2  •  Pj  is  initial  pressure  (torr),  7  is  j^sec, 
laboratory  time  d  is  100  times  initial  fraction  of  CO. 
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Tb~  same  data  as  in  (a)  plotted  as  function  of  rpj^  (  a  very 
similar  scatter  is  obtained  if  rp^  is  used  as  ordinate). 
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Fig.  4c  Incubation  times  from  ultraviolet  experiments  at  2478^.  The 
solid  line  is  taken  from  Fig.  4a. 
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By  analogy  with  the  chemistry  of  air  the  following  reactions 
require  consideration: 


CO  +  M 

^  C  +  O  +  M 

(1) 

C  +  CO 

-  c2  +  o 

(2) 

O  +  CO 

~  °2  +  C 

(3) 

C2  +  M 

~  C  +  C  +  M 

(4) 

°2  +  M 

~  O  +  O  +  M 

(5) 

Of  these,  only  the  rate  of  (5)  is  known.  The  list  is  not  exhaustive  and 

other  reactions  are  possible.  At  the  high  temperatures  of  these  experiments, 

dissociation  reactions  predominate  and  it  is  felt  unlikely  that  triatomic 

species  will  be  important  routes  in  dissociation.  For  instance,  the  reaction 
_  9 

CO  +  O^  •*-  CO^  +  O  has  been  measured  but  in  order  to  become 
atomized  any  CO^  will  have  to  dissociate  again.  This  matter  of  triatomic 
molecule  formation  will  be  more  fully  discussed  later  on. 

The  correlations  of  incubation  time  shown  in  Fig.  4  are  two  of  those 

tried.  A  plot  of  log  Tp^  vs  l/T  and  Tp$  vs  l/T  show  similar  scatter  and 

1  /2  . 

it  was  found  that  log  T  p<p  gives  the  best  correlation.  The  temperatures 
used  are  those  assuming  frozen  chemistry,  but  full  vibrational  and  rotation¬ 
al  relaxation.  In  the  dilute  mixtures  used  the  density  jump  at  the  shock 
front  is  always  close  to  4  and  the  parameter  p^</>  is  a  measure  of  the 
CO  concentration. 

The  incubation  results  may  be  explained  on  chemical  grounds,  but 
it  is  necessary  to  first  rule  out  shock  attenuation  as  the  cause.  Attenuation 
will  have  two  main  effects  on  the  observations.  The  increasing  temperature 


behind  the  front  will  cause  more  rapid  reaction,  particularly  in  view  of  the 

very  high  activation  energy  of  dissociation  of  CO.  The  boundary  layer 

build-up  will  also,  even  in  the  absence  of  a  temperature  increase,  by 

shortening  the  time  scale  in  laboratory  co-ordinates  cause  an  apparent 

increase  in  reaction  rate.  But  the  small  fraction  of  CO  in  the  mixture  will 

1  /2 

not  alter  the  aerodynamics  much  and  the  scaling  found  depending  on^  ' 
is  not  to  be  expected  unless  it  is  due  to  the  chemistry. 

Reactions  of  the  type  O  +  CO  -*  O  +  C  +  O  and  C  +  CO  —  C  +  C  +  O 
will  have  an  auto  catalytic  effect  but  with  concentrations  of  1%  or  less  of 
CO  integration  of  these  equations  together  with  (1)  showed  that  O  or  C  will 
need  to  be  about  a  factor  of  1000  faster  than  argon  in  the  dissociation  of 
CO  to  give  the  observed  incubation  behavior  even  approximately.  This 
factor  is  unlikely;  comparison  with  oxygen  and  nitrogen  dissociation  shows 
the  effect  to  be  about  a  factor  of  25  or  35. 

It  thus  appears  that  some  diatomic  intermediate  is  formed  and  this 

causes  the  turn  over  from  an  initially  slow  rate  to  a  much  faster  one. 

* 

There  are  three  a  priori  possibilities,  O^,  C^  and  CO  where  the  asterisk 
refers  to  any  or  all  of  the  triplet  states  of  CO  which  dissociate  to  normal 
products.  ^  In  order  to  affect  the  dissociation  rate  to  give  an  incubation 
effect  as  observed,  it  is  necessary  that  the  intermediate  should  have  a 
higher  cross  section  for  dissociation  than  the  original  CO.  The  Keck- 
Carrier^’  ^  theory  of  dissociation  may  be  used  to  estimate  that  the  recom¬ 
bination  rate,  kr,  for  (1)  with  CO  in  the  ground  state  is~10  ^  cm^  sec  * 

3  -34  #  -33 

while  for  C^  in  the  X  FI  state  it  is  «  7x10  and  for  CO  it  is«  10 
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The  rate  is  experimentally  determined  to  be  close  to  10  .  Because 

of  the  low  equilibrium  concentration  of  in  dissociating  CO  mixtures 
unless  raction  (3)  is  made  very  fast  (»gas  kinetic  cross  sections)  cannot 
be  the  intermediate. 

Triplet  States  of  CO 

The  observation  of  the  Cameron  band  (0,  1)  at  2155$  indicates  that 

3  1 

triplet  states  of  CO  are  present,  the  transition  involved  being  a  FI  —  X  2. 

Davies*^  has  reported  the  observation  of  the  triplet  bands  at  6434$  , 

3  3 

due  to  d  A  —an,  and  3A  bands  in  the  ultraviolet.  Although  they  were 
carefully  sought,  these  triplet  and  3A  bands  were  not  positively  identified 
in  these  experiments  by  photographic  means.  A  series  of  experiments  was 
run  using  photoelectric  detection.  The  CO  triplet  radiation  will  occur  be¬ 
tween  4000$  and  9000$  ,  but  the  spectrum  is  dominated  between  4000  and 
6000$  by  C2  Swan,  and  CN  red  will  occur  between  8000  and  9000$  . 
Observations  were  made  from  6000  -  12,000$  using  Corning  filters  and  a 
SI  surface  photomultiplier.  Several  combinations  were  used.  Some  were 
broad  band  and  some  narrow  (500$  ).  Because  of  C^  and  CN  and  continuum 
radiation,  it  was  not  possible  to  positively  identify  any  radiation  as  due  to 
the  CO  triplet  states  in  this  region. 

3 

The  mode  of  excitation  to  the  a  FI  state  is  a  matter  of  importance 

1  3 

but  little  knowledge.  Estimates  have  been  made  that  the  Cameron  bands 

3  11 

are  around  a  factor  of  10  weaker  than  the  fourth  positive,  A  n  —  X  2. 

Recent  work^  has  shown  the  lifetime  of  the  A^n  state  to  be  around  10  ® 

3  -5 

sec  and  so  the  radiative  lifetime  of  the  a  FI  is  probably  about  10  sec. 

*A11  rates  in  this  paper  will  be  given  in  particles  per  cubic  centimeter  per 
second  units . 
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Using  Hansche's  observations  on  the  quenching  of  the  Cameron  bands  in  a 
discharge  it  may  be  concluded  that  the  cross  section  for  the  process 
CO(*2 )  +  CO(3n)  -®  CO(^2 )  +  CO(^£ )  is  close  to  gas  kinetic. 

It  is  unlikely  that  argon  will  have  a  substantial  cross  section  for  the 
3  -2 

deactivation  of  CO(a  IT),  it  is  probably  less  than  10  of  gas  kinetic  and 
thus  unimportant  in  1%  CO  mixtures,  compared  with  CO  itself.  Atoms, 
either  C  or  O  are  not  expected  to  have  cross  sections  significantly  greater 
than  gas  kinetic,  thus  in  a  1%  mixture  where  less  than  10%  of  the  CO  has 
dissociated  atoms  will  still  be  less  important  that  CO(*£)  in  the  deactiva¬ 
tion.  Radiative  deactivation  will  always  set  an  upper  time  limit  on  the 

3 

achievement  of  a  steady  state  of  CO(a  II)  if  the  major  process  involves 

-5 

CO  molecules  and  not  atoms.  This  time  is  of  the  order  10  sec  or  about 
3  fisec  in  laboratory  coordinates.  The  incubation  times  were  generally 
much  longer  than  this,  indicating  that  there  is  no  connection  between  the 
incubation  time  and  CO  if  Refs.  (13)  and  (15)  are  correct  in  their  findings. 

The  incubation  time  appears  most  likely  to  be  associated  with  the 
formation  and  dissoc:ation  of  C^.  It  was  remarked  earlier  that  the  C^  con¬ 
centration  reached  a  peak  after  the  end  of  the  incubation  time.  The  absolute 
concentration  of  C^  is  always  much  less  than  that  of  CO  and  so  at  all  times 
the  C  and  O  atom  concentration  are  nearly  equal. 

The  end  of  the  incubation  time  is  noted  experimentally  to  occur  when 
[  C 2]  T  88  ^2  [CO]/l0  where  is  the  equilibrium  constant  of  reaction  (2) 
and  the  time  required  to  achieve  this  is  t  w[C_]  /d[C~]/dt.  From  reactions 

(2)  and  (4)  d[  C2]/dt  =  k2[CO]  [C  ]  -  k2  [C2]  [O  ]  -  k4M  [C2  ]  ,  where  k  and  k" 
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refer  to  the  forward  and  backward  rate  constants,  and  at  early  times  the 
second  term  is  not  important.  The  first  and  third  terms  are  nearly  equal 
and  the  difference  is  probably  of  order  [CO]  [C]/lO.  The  average  of  [C] 
at  t  is  k^CO]  [M]t/2.  Thus 

7  k_k.  [CO]  2[M]  t  dt  k,k.  [CO]2  [M]r2 

[C2b  ■  J  10  *"'2 -  »  1  40 - 


hence 


4K. 


7 


k2kj  [M]  [CO] 


since  [CO  ]  is  taken  as  constant  during  the  incubation  ffme. 

An  expression  of  similar  form  but  different  value  is  obtained  by 
noting  that  the  end  of  the  incubation  time  will  be  apparent  when 


d[C2]/dt  *  2k^  [CO]  [M  ] 

from  which 

lc2]T  [M]k4  *  2k,[COJ|M], 

as  before 

[C 2]T  *  kik2tC°]2  T2/40, 

hence 

2  80 
7  «  - 

k2k4  [M][CO] 

These  two  expressions  for  the  incubation  time  are,  of  course,  related  by 
the  ratio  of  k^  to  k^.  [M  ]  is  mainly  argon  and  is  closely  proportional  to 
PjA,  while  [CO]  is  proportional  to  p^A.  A,  the  density  jump  across  the 
shock,  does  not  very  greatly  and  so  7  p^  0  »  K2/k2kj.  It  is  also 
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necessary,  as  mentioned  before,  for  the  production  of  carbon  atoms  via 

reactions  (2)  and  (4)  to  exceed  the  production  from  reaction  (1)  in  order  to 

observe  the  change  of  slope  in  the  CO  trace.  Thus,  if  reaction  (2)  is  fast 

enough  to  keep  in  local  equilibrium,  the  condition  for  a  clearly  observed 

change  in  slope  is  [CO]  [M  ]  >  [CO]  [M  ]  kj. 

2  2 

In  Fig.  4a  is  plotted  1/2  log  t  Pj  <j>  vs  1 /T.  The  slope  of  this  line 

1  /2 

should  be  given  by  (K^/k^k^)  .  has  an  activation  energy  of  113  K 

cal/mole  while  k^  will  probably  have  a  similar  activation  energy,  hence  the 

1  /2 

temperature  variation  of  the  incubation  time  will  be  close  to  that  of  (l/k^)  ' 

or  128  K  cal/mole.  The  value  obtained  from  the  data  is  135  +  15  K  cal/mole 

which  is  considered  a  satisfactory  agreement. 

POSSIBLE  ALTERNATE  INCUBATION  MECHANISM 

In  a  previous  communication,  ^  the  author  argued  that  the  incubation 

phenomenon  was  more  probably  to  be  associated  with  CO  before  becoming 

1  5 

aware  of  the  implications  of  reference.  On  the  basis  of  the  data  presented 

in  this  paper  it  is  not  in  fact  possible  to  separate  out  the  possible  effect  of 
$ 

CO  if  it  should  turn  out  to  have  a  much  larger  lifetime  than  that  assumed 

-5  1  3 

(10  s)  or  if  CO(  2)  has  a  very  much  lower  cross  section  than  0(  P)  or 

C(^P)  for  the  excitation  to  CO  . 

Extraneous  Effects  on  the  Observations 

At  the  high  temperatures  used  in  these  experiments  there  are 
several  effects  which  can  interfere  both  with  the  reaction  mechanism  and 
with  the  observations.  Fig.  3e  demonstrates  the  effect  of  ionization  of 
carbon  on  the  infrared  signal.  The  CO  emission  will  tend  to  zero  as  the 
CO  dissociates,  but  the  continuum  radiation  (free  bound  plus  free  free)  due 
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to  ion-electron  and  neutral-electron  interaction  will  rise  as  the  gas  ionizes 
at  these  high  ,emperatures.  Because  of  this  effect,  it  is  not  possible  to 
follow  the  CO  concentration  to  the  completion  of  the  reaction,  and  even  the 
initial  portions  of  the  trace  are  distorted. 

Another  effect  which  might  interfere  with  a  straightforward  analysis 

of  the  data  could  be  the  presence  of  argon  metastable  'oecies.  These  have 

17 

been  shown  by  Petschek  and  Byron  to  be  important  in  the  initial  ionization 
of  argon  by  atoms.  They  occur  at  an  excitation  energy  of  11.5  ev,  which  is 
fairly  closely  resonant  with  the  dissociation  energy  of  CO,  11.1  ev.  To 
check  this  possibility  some  runs  were  made  with  neon  as  the  working  gas. 

The  lowest  metastable  levels  of  neon  lie  at  16.  6  ev,  and  accidental  inter¬ 
ference  with  the  dissociation  will  not  occur.  The  records  show  that  there  is 
little  difference  between  argon  and  neon  as  the  collision  partner.  Because 
the  neon  has  to  be  driven  by  combustion,  attenuation  of  the  shock  wave 
gave  rather  larger  scatter  in  the  measured  incubation  time. 

Effects  of  Impurities  on  the  Incubation  Time 

The  effect  of  impurities  on  the  incubation  time  was  studied  in  a 
series  of  experiments.  The  additives  tested  were  oxygen,  added  to  replace 
5,  20  and  50%  of  the  original  CO,  nitrogen  and  carbon  tetrachloride  to  replace 
5%  of  the  CO.  Carbon  per  se  cannot  be  added  to  the  CO  and  the  addition  of 
carbon  tetrachloride  was  considered  as  the  easiest  means  of  adding  carbon 
atoms.  Unfortunately,  it  will  give  four  times  as  many  chlorine  atoms  and 
while  it  might  be  supposed  that  chlorine  atoms  will  act  as  an  inert  diluent, 
a  series  of  experiments  should  be  run  to  check  this  point.  Under  the 
conditions  of  these  experiments  both  O^  and  CCl^  will  be  dissociated  in 
much  less  than  the  incubation  time. 
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The  results  are  shown  in  Fig.  5.  The  three  solid  lines  are  taken  from 
Fig.  4a  and  are  for  reference  only.  It  will  be  noted  that  there  is  a  larger 
scatter  and  that  nitrogen  has  no  significant  effect.  There  appears  to  be  a 
general  slight  shortening  of  the  incubation  time  with  carbon  tetrachloride 
and  a  lengthening  with  oxygen  additions.  While  the  results  are  not  con¬ 
clusive,  they  suggest  that  carbon  atoms  do  shorten  the  incubation  time. 

While  the  effect  of  upon  the  incubation  time  was  being  investigated, 
its  effect  on  the  CN  radiation  was  also  studied.  As  might  be  expected, 
addition  of  to  the  mixture  strongly  enhanced  the  CN  emission  but  did  not 
change  any  other  radiation. 

The  Initial  Reaction  Rate 

The  majority  of  the  argon  diluted  experiments  did  not  show  an  initial 

slope  measurably  different  from  zero.  On  the  other  hand  runs  in  neon 

mixtures  tended  to  give  a  finite  initial  slope.  This  difference  in  behavior 

indicates  strongly  that  even  at  early  times  ionization  is  occurring  and  a 

rising  component  of  ionization  continuum  is  superposed  on  the  falling  CO 

signal.  The  substitution  of  neon  for  argon  will  at  once  tend  to  give  lower 

ionization  levels  and  since  it  also  has  a  much  smaller  cross  section  for 

1 8 

neutral  Bremstrahlung  the  interference  with  the  CO  signal  will  be  reduced. 

Because  of  the  auto  catalytic  effect,  it  is  to  be  expected  that  the 
initial  slope  showed  gradual  increase  and  taking  into  account  the  noise  on  the 
traces  and  the  fact  that  aerodynamic  effects  on  the  temperature  can  drastically 
change  the  reaction  rate  it  is  not  expected  that  the  initial  rate  data  which  was 
obtained  and  is  presented  in  Fig.  6  is  particularly  reliable.  In  fact,  it  is 
probably  not  better  than  about  a  factor  of  three.  The  scaling  found  to  give 
the  best  representation  of  the  data  was  d  [CO]/dt  =  -  kj  [M]  [CO], 
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Fig.  5  Effects  of  impurity  addition  on  incubation  time.  The  lines  are 
taken  from  Fig.  4a. 
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Fig.  6  CO  dissociation  rates  obtained  from  initial  slopes.  The  line 
marked  D  is  Davies'  ratelZ  and  PCW  is  from  Presely, 
Chackerian  and  Watson.  19 
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The  line  which  was  fitted  to  the  data  by  eye  and  anticipated  the  low 

temperature  observations  would  give  rates  which  were  to  fast,  is  represented 

by  kj  =  1  x  lO'10^)1/2  exp(-D/RT)  cm3  s  *  which  is  equivalent  at  8000°K  to  a 

recombination  rate  constant  of  4  x  10  cm^  s  ^ .  Also  plotted  in  Fig.  6 

1 2 

are  lines  due  to  Davies  for  Ar  plus  CO  mixtures  and  Presely,  Chackerian 
19 

and  Watson  for  pure  CO.  From  the  way  in  which  these  authors  obtained 
their  rates  they  should  be  more  properly  compared  with  what  will  be  called 
the  second  reaction  rate  in  this  paper. 

Some  of  the  data  presented  by  Presley,  Chackerian  and  Watson  for 
the  infrared  emission  from  dissociating  pure  CO  indicate  the  presence  of 
an  incubation  time  although  the  authors  do  not  comment  on  this.  Their 
shock  front  temperatures  were  around  20,  000°K.  When  plotted  on  Fig.  4a, 
extrapolating  the  solid  line  given,  it  is  found  that  the  CO- CO  incubation  time 
is  about  50%  longer  than  the  CO-Ar  (extrapolated)  value.  This  difference 
lies  within  the  experimental  error  particularly  since  shock  curvature 
introduces  a  significant  uncertainty  at  the  low  pressures  used  in  the  pure 
CO  work.  It  is  improbable  that  CO-CO  collisions  are  less  effective  than 
CO-Ar  and  the  data  of  Presley,  Chackerian  and  Watson  is  interpreted  here 
to  mean  that  CO-CO  and  CO-Ar  have  comparable  cross  sections  in  the 
dissociation  process. 

Second  Reaction  Rate 

A  temperature  change  of  100°K  can  cause  a  change  in  reaction  rate 
of  20%  in  the  range  of  interest.  Such  a  temperature  drop  can  be  caused  by 
a  fairly  small  degree  of  reaction.  Acting  in  the  opposite  direction  attenuation 
and  boundary  layer  growth  will  actually  cause  a  temperature  rise  which  may 
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more  than  off  set  the  drop  due  to  reaction  and  given  an  apparently  faster 
reaction.  The  data  presented  in  Fig.  7  are  affected  by  these  considerations, 
in  particular  those  taken  at  the  low  temperature  end.  In  these  the  slope  was 
measured  some  distance  behind  the  shock  front  and  from  the  known  attenuation 
alone  the  observed  rates  at  6500°K  are  probably  a  factor  of  two  high.  At 
the  high  temperature  end  of  the  range,  ionization  begins  to  cause  lower 
rates  to  be  observed  and  thus  the  trend  in  the  data  will  be  to  give  a  lower 
activation  energy  than  the  true  value. 

It  will  be  noted  that  the  infrared  data  acquired  in  the  present  experi¬ 
ments  is  in  quite  good  agreement  with  that  of  Davies,  although  the  inter¬ 
pretation  is  entirely  different. 

In  view  of  the  arguments  developed  concerning  the  incubation  time, 
it  is  clear  that  this  second  reaction  rate  probably  does  not  have  a  large  com¬ 
ponent  of  CO  dissociation,  rather  it  is  determined  by  reactions  (2),  (3),  (4) 
and  (5)  and  while  it  can  be  fitted  with  a  kinetic  rate  expression  referring  to 
[  M  ]  and  [  CO  ] ,  such  procedure  has  only  empirical  value. 

THE  RATIO  OF  THE  INITIAL  TO  SECOND  REACTION  RATE 

Experimentally  it  is  found  that  there  is  about  a  factor  of  6  between 
the  first  and  second  rates  at  8000°K.  It  is  also  estimated  that  the  C^  con¬ 
centration  during  the  final  stages  of  the  reaction  is  about  a  factor  of  two  be¬ 
low  its  equilibrium  value.  The  rate  of  disappearance  of  CO  via  reaction  (1) 
is  d  [C^]  dt  =  [CO]  =  -  k^  [CO]  [M]  while  via  reactions  (2)  and  (4)  the  rate 

is  given  by  [COJ2  =  -  k2[CO]  [C]  +  k2  [C2][0]  and  [<~2]  =  -  [CO]2  -  [C2] 

•  •  • 

[M]  but  since  [CO]  »[C2],  [CO]2  ~  -  [C2]  [M]  k^  and  substituting  [C2]  85 
K2[CO]/2,  (the  exp  erimentally  observed  limit)  the  result  is  [CO]2  ~  K2 
[CO]  [M]  k^/2.  Thus  by  rearranging  k^  ~  10  k^K^/k^K^  and  at  8000°K, 
k'  *  10  k" 
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CO  dissociation  rates  obtained  from  the  final  slopes 
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liehavior 

Of  the  possible  products  of  shuffle  reactions  only  was  observed. 

1  3 

In  the  photographic  spectra  taken,  both  2  and  n  were  identified.  The 

3  3 

photoelectric  studies  concentrated  entirely  on  the  Swan  bands  (  IT  — *  FI). 

While  in  reaction  (2)  it  is  reasonable  to  suppose  the  C  and  O  atoms  to  be  in 
3  1 

their  P  ground  states,  and  CO  in  the  X  2  state,  it  is  not  obvious  what  state 

1  3 

the  will  be  in.  The  X  2  and  X  n  states  are  very  close  in  energy  and  in 

the  data  reduction  it  has  been  assumed  that  full  equilibrium  exists  between 

these  two  states.  It  was  also  assumed  that  the  observed  emission  was  in 

local  equilibrium  and  absolute  concentrations  of  were  obtained  by  the  use 

20 

of  a  previously  obtained  f  number.  In  some  of  the  experiments  at  low  CO 
concentrations  and  low  initial  pressure,  it  is  possible  that  the  radiation  will 
be  collision  limited  and  so  the  concentration  deduced  from  emission  will  be 
lower  than  the  true  value.  At  most  this  effect  will  not  exceed  a  factor  of  two. 

The  emission  in  the  Swan  bands  showed  a  gradually  increasing  slope 
forming  a  toe  at  the  shock  front.  Within  a  few  microseconds  the  rise  be¬ 
came  effectively  linear.  The  radiation  went  through  a  peak,  which  occurred 
after  the  end  of  the  incubation  time  and  then  declined  with  the  same  profile 
as  the  CO  indicating  that  the  C^  and  CO  were  in  steady  state.  Towards  the 
end  of  the  dissociation  the  C^  is  a  more  useful  observable  than  the  infrared 
which  becomes  dominated  by  Bremstrahlung. 

The  linear  part  of  the  C^  rise  was  measured  for  a  range  of  conditions 
and  the  results  given  in  Fig.  8.  Using  the  mechanism  of  reactions  (1),  (2) 

and  (4),  putting  C^  in  steady  state  and  neglecting  back  reactions  it  may  be 

2 

estimated  that  d  [C 2  ]/dt  «2k^k2  [CO]  /k^.  The  use  of  any  other  power  of 
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Fig.  8  Cz  formation  rate  during  incubation  time.  The  ordinate  is  in 
mV/jisec,  C  is  a  calibration  factor,  G  is  the  solid  angle 
subtended  by  the  gas  at  the  spectrometer,  is  the  fraction 

of  C 2  molecules  radiating  in  the  band  pass.  By  the  use  of  the 
f  number  in  Ref.  20,  d  [C2  ]  / dt  (particles/cc  sec)  is  given 
by  3.  15  x  1017  x  ordinate  for  p ^  =  1  torr  ,  $  -  1.0.  Figures  9 
through  16  present  comparisons  between  calculations  and 
experiment.  All  times  are  given  in  laboratory  coordinates. 
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The  activation  energy  found 


[  CO  jgave  much  greater  scatter  of  the  points. 

200  k  cal/mole,  is  consistent  with  the  ratio  k^k^/k^.  The  reaction 

co  +  co  -  c2  +  o2  (6) 

has  an  endothermicity  of  250  K  cal/moie  and  if  an  estimate  of  0.  28 
21 

(D  +  T>„  )  is  added  to  this,  a  total  activation  energy  of  about  320 
°2  C2 

K  cal/mole  is  expected.  By  the  use  of  pre-exponential  factors  in  the  rate 
expression,  it  is  possible  to  fit  the  slope  with  the  form  k  <x  (E/RT)n  exp 
(-E/RT)  with  E  =  320  K  cal/mole  and  n  =  9.  For  a  collision  of  two  CO 
molecules  the  greatest  value  of  n  which  is  classically  allowed  is  4  and  so 
the  above  reaction  cannot  contribute  to  the  C2  formation  in  a  major  way. 

The  other  shuffle  reaction,  forming  02  via  reaction  (3)  was  not 
directly  observed.  Its  effect  was  deduced  from  some  experiments  in  which 
02  was  added  to  the  CO  mixture.  With  equal  amounts  of  02  and  CO  diluted 
with  98%  Ar,  there  was  an  initial  reaction  rate  almost  obliterating  the 
incubation  time  since  there  was  no  point  at  which  the  slope  obviously  changed. 
Fig.  1  6a  .  These  experiments  will  be  more  fully  discussed  in  the  section  on 
computation. 

Other  Reactions 

It  has  already  been  mentioned  that  many  reactions  may  be  written 
down  for  this  system.  For  instance 

CO  +  CO  —  CzO  +  O  (7) 
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(C^O  has  recently  been  discovered  as  a  result  of  the  action  of  C  atoms  on 
22 

CO  in  a  matrix)  and  its  analogue 

CO  +  CO  —  C02  +  C  (8) 

are  also  candidates  for  the  removal  of  CO.  If  the  latter  reaction  should 

be  the  producer  of  carbon  atoms  required  for  reaction  (2)  rather  than 

reaction  (1),  then  the  incubation  time  parameter  should  be  t  p  </>  rather  than 

1  /2 

T  p  <j>  '  and  the  activation  energy  of  the  plots  of  Fig.  4  should  be  around 
65  K  cal/mole,  neither  of  which  is  in  accord  with  the  experiments.  These 
considerations  indicate  that  the  above  reactions  are  not  dominant.  Such 
is  the  scatter  of  data  that  it  cannot  be  shown  whether  or  not  such  reactions 
do  occur  to  a  minor  extent. 

C^O  may  also  participate  in  the  reaction  scheme  by  producing  the 
C^  by  a  reaction  such  as 


c2o  +  C  —C2  +  CO 


(9) 


and  if  the  main  removal  of  C20  is  the  decomposition 


C20  +  M—CO  +  C+  M 


(10) 


The  heat  of  dissociation  of  reaction  (10)  is  not  known  for  certain,  but  has 
been  estimated  as  44  K  cal/mole.  This  value  is  probably  low  but  not  of 
great  consequence  for  the  following  analysis.  The  steady  state  analysis 
already  applied  to  the  formation  of  C2  from  the  reactions  (1),  (2)  and  (4) 
gives 

k  [CO][C]  .  2k,k  [CO]2 

[C2J«  — -  and  [C  ]  «  — - 

MM]  k4 
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a  similar  analysis, 
( 7)  and  (10)  gives 


putting  both  [C^]  andfC^O]  in  steady  state,  using  reaction 


[c2o] 


k?  [CO]2 
k1()[M] 


and 

k  [C  0][C]  k_  [CO]  ^  [C  ] 

[C  ]  =  -5 - ? -  .  thus[c  ]  =  _Z - 

k4[M]  k4k1Q[M]^ 

and  if  the  C  atoms  are  produced  mainly  by  reaction  (8)  we  obtain 

k_ka  [co]4 

[C2]  A  =  — - T 

k4k1Q[M]2 

and  if  (1)  is  the  reaction  to  produce  carbon  atoms, 

kjk7  [CO]3 

k4k10 

Neither  of  these  expressions  fits  the  data  of  Fig.  8,  either  in  the  concentration 
or  slope  dependences,  giving  respectively  56  and  181  K  cal/mole  for  the 
activation  energy.  By  adjusting  the  heat  of  reaction  of  10  to  84  K  cal/mole, 
a  value  which  is  probably  too  high,  the  slope  may  be  made  to  fit,  but  the 
concentration  dependence  remains  at  vJ.riance  with  the  data. 

For  these  reasons,  it  is  not  felt  that  reactions  (7)  through  (10)  have 
any  importance  in  the  high  temperature  dissociation  of  carbon  monoxide. 
Because  of  the  unimportance  of  in  the  dissociation  of  CO,  the  possible 
triatomic  reactions  leading  to  it  may  be  ignored. 

Impurity  Effects 

It  has  been  commented  earlier  that  nitrogen  impurity  in  the  gas  gave 
rise  to  CN  radiation.  It  is  probable  that  CN  too  will  be  formed  in  shuffle 
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reactions  and  might  thus  act  as  another  intermediate  in  the  CO  dissociation 
path.  Nitrogen  impurity  added  in  amounts  up  to  5%  of  the  carbon  monoxide 
content  has  no  observable  effect  on  the  second  rate,  neither  did  a  similar 
addition  of  oxygen  or  of  carbon  tetrachloride.  The  effect  of  impurity 
addition  on  the  initial  slope  was  also  not  significant  in  the  case  of  nitrogen 
and  oxygen  but  the  carbon  tetrachloride  experiments  did  show  about  a  factor 
of  two  enhancement  of  the  rate  at  temperatures  below  8000°K.  There  is 
very  little  data  above  that  temperature. 

COMPUTER  STUDIES 

Of  the  reactions  mentioned  as  probably  important,  i.  e.  ,  (1)  through 
(5),  only  one,  the  oxygen  dissociation  rate,  has  been  measured.  All  of  the 
others  are  unknown.  Parametric  computer  studies  were  undertaken  to  find 
a  plausible,  self-consistent  set  of  rate  constants  which  will  predict  both 
the  CO  and  C^  histories  foi  a  wide  range  of  conditions.  The  experiments 
all  lie  in  the  range  of  temperature  between  about  6000°K  below  which  they 
are  too  slow  and  shock  attenuation  effects  become  troublesome  and  10,  000°K 
where  the  experimental  time  is  very  short  and  ionization  gives  spurious 
signals. 

Because  of  the  complexity  of  the  reaction,  it  was  not  considered 
meaningful  to  try  to  establish  the  temperature  variation  of  the  backward 
rate  constants  (all  rates  will  be  given  in  the  backward  direction  as  written* 
i.  e.  ,  in  the  exothermic  direction,  by  relating  the  forward  and  backward 
rates  through  the  equilibrium  constant).  Rather,  the  attempt  will  be  made 
to  establish  the  rates  to  a  close  order  of  magnitude  in  this  temperature 
range. 
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Two  different  programs  have  been  used.  One  is  based  on  the  steady 
state  approximation  for  the  intermediate  species  C_,  O  ,  the  rate  of  change 

w  £ 

of  whose  concentration  is  very  much  less  than  that  of  C,  O,  or  CO. 

The  final  equilibrium  conditions  were  found  from  an  equilibrium 
program  due  to  Los  Alamos  Laboratory,  neglecting  ionization,  thus  the 
final  degree  of  dissociation  was  found.  The  reaction  was  divided  into  100 
equal  steps  of  dissociation  and  the  average  concentration  of  each  of  the 
species  found  at  each  step.  An  estimated  set  of  rate  constants  was  used 
and  thus  the  time  required  for  each  given  change  in  a  ,  the  degree  of 
dissociation  was  found.  A  running  integral  of  the  time  was  maintained,  con¬ 
verted  to  lab  time  via  the  relation  At^  =  At  /p~  ,  where  p  is  the  average 
density  ratio  across  the  shock  at  point  of  interest. 

Boundary  Layer 


Many  of  the  experiments  were  followed  well  into  the  test  gas 
(laboratory  times  >20  p,  s).  In  these  experiments  boundary  layer  growth 
is  important.  This  leads  to  four  main  effects.  The  boundary  layer  reduces 
the  path  length  through  the  hot  gas  and  sheaths  this  hot  core  with  cold  gas. 
For  the  CO  infrared  and  C^  visible  radiation,  this  will  not  affect  the 
radiation  as  seriously  as  the  other  concomitants  of  the  boundary  layer  which 
comprise  a  decrease  in  particle  velocity  with  respect  to  the  shock  front, 
leading  to  an  increase  of  enthalpy  of  a  particle  behind  the  shock  and  to 
an  increase  in  particle  time  which  is  now  given  by  t  =  t.  p  F  (t)  where  F  (t) 

O  » 

is  greater  than  unity.  These  effects  will,  in  the  absence  of  corrections, 
give  rise  to  deduced  reaction  rates  which  are  too  great.  Towards  the  end 
of  the  test  gas  the  effect  in  time  alone  can  be  a  factor  of  about  5.  Pari 
passu  with  the  boundary  layer  growth  a  deceleration  of  the  shock  occurs. 
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This  also  leads  to  progressively  higher  enthalpy  behind  the  shock.  Mirels  has 
given  expressions  applicable  to  the  center  stream  line  for  the  time  and 
particle  velocity  transformation. 


which  gives 


U 

U 


e 

eo 


1  - 


2PZ/P 1 


£  m 


=  £  n  (1 


1/2 


£  rrv 


+  {T~cn] 


1/2 


where  is  particle  velocity  with  respect  to  the  shock,  U0O  is  the  velocity 
at  the  shock  front.  £  is  the  distance  behind  the  shock  front,  £m  is  the 
maximum  length  of  the  test  slug  and  is  obtained  from  Mirels'  curves 
assuming  a  laminar  boundary  layer.  xg  =  t^Us  and  £  =  t^Ug  where  Ug  is  the 
shock  velocity.  The  shock  attenuation  effect  on  enthalpy  was  estimated  by 
assuming  a  linear  enthalpy  rise  behind  the  shock  with  distance.  The  shock 
velocity  corresponding  with  the  gas  at  the  contact  face  was  estimated  from 
the  velocity  measurements  made  along  the  tube.  We  thus  obtain 


AHe  =  l/2  M(Ug^  -  Ue^)  ,  M  is  the  molecular  weight 
U  =  U  +  £/£  (U  -  U  )  ,  x  denotes  the  contact  surface. 

G  S  XX  S 

In  reality,  the  situation  is  more  complicated,  since  wave  interactions  are 

involved  and  not  just  enthalpy  addition.  These  effects  have  never  been 

properly  treated  and  lie  beyond  the  scope  of  these  studies. 

3 

The  Avco  stream  tube  program  has  also  been  used  to  study  the 
reaction  system,  using  as  input  the  initial  estimates  obtained  from  the 
steady  state  program. 


-33- 


As  is  to  be  expected,  the  effect  of  the  boundary  layer  increases 
progressively  with  distance  behind  the  shock  and  about  halfway  through 
the  test  gas  the  correction  is  about  a  factor  of  two.  Since  the  cor¬ 
rection  is  somewhat  uncertain  itself,  there  is  little  value  in  analyzing  the 
data  beyond  this  point.  The  effect  of  the  boundary  layer  on  the  incubation  is 
within  the  experimt  .ital  error  of  measuring  the  time. 

Reaction  Rates  (Estimates  from  Data) 

Since  only  the  rate  of  reaction  (5)  is  known,  it  is  required  to  find  the 
rates  of  the  other  four.  In  order  to  do  this  it  is  necessary  to  make  estimates 
of  what  is  likely  to  be  the  upper  bounds  of  the  rates  so  that  the  computations 
are  not  forced  to  fit  the  data  by  a  set  of  unrealistic  assumptions.  As  a 
starting  point,  the  estimates  made  from  the  experiments  were  used.  It  has 
already  been  shown  that  k^  «  10  k^  and  in  conjunction  with  the  rise 

*  *5  1C 

experiments,  from  which  [C^]  /  [CO  ]  =  1.9  x  10  at  8000°K,  k^  may 

be  estimated,  since  1.9  x  10  =  2k^  k^  /k^  from  the  expression 

given  earlier.  Thus  k^  ~  3  x  10”*®. 

-  34 

The  initial  slope  measurements  gave  k^  «  4x  10  and  the  second 

33  o 

slope  measurements  of  the  CO  traces  yield  k^  «  4  x  10  at  8000  K 
neglecting  the  contribution  of  k^  and  taking  [C  a  l/2  K 2  [CO  ]. 

Because  of  the  approximations  involved  in  deriving  the  various 
kinetic  expressions,  it  is  not  expected  that  the  rates  obtained  from 
them  are  accurate.  The  initial  input  to  the  calculations  was  taken  as 
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kj  =  10  =  3  x  10  k^  =  3  x  10  ^  and  thereafter  these  rates  were 

adjusted  to  obtain  as  good  fits  as  possible  over  a  range  of  experiments. 
Samples  of  data,  along  with  the  computer  analysis,  are  shown  in 
Figs.  9-16;  these  are  discussed  in  appropriate  places  below,  but  details 
are  given  in  the  figure  captions. 

By  assuming  arbitrary  values  for  k^  it  was  found  that  reaction  (3) 

-  8 

had  no  perceptible  effect  unless  >  10  .  This  is  an  unreasonably  high 

value  and  so  the  participation  of  as  an  intermediate  was  ignored  in  CO 
and  noble  gas  experiments.  In  this  determination  attention  was  paid  to 
both  the  and  CO  profiles.  One  of  the  main  difficulties  was  experienced 
in  the  matching  of  the  CO  profiles.  It  has  already  been  noted  the 
Brems strahlung  probably  obscures  the  initial  decay  of  the  CO  and  affects 
the  observed  later  decay  as  well.  The  rising  tail  of  the  infrared  signal  in  a 
run  such  as  shown  in  Fig.  10a  was  rather  arbitrarily  extrapolated  linearly 
back  towards  the  shock  front  and  the  CO  signal  taken  as  the  difference  be¬ 
tween  this  line  rather  than  the  base  line  and  the  trace.  Figure  10b  also  shows 
the  CO  data  reduced  with  respect  to  the  base  line.  The  difference  is  greater 
than  the  uncertainty  of  reading  the  trace  which  is  also  indicated  in  Fig.  9(d). 

It  is  clear  that  in  future  work  the  Brems  strahlung  should  be  measured  in¬ 
dependently  in  order  to  allow  for  its  effect  on  the  CO  channel. 

While  unimportant  in  CO  and  noble  gas  mixtures,  reaction  (3)  is  of 
importance  when  much  oxygen  is  present. 

It  will  be  noted  that  there  are  two  major  discrepancies  in  the  calculated 
and  observed  C ^  histories.  In  the  10%  CO  mixtures  the  C^  level  computed 
was  higher  than  that  observed  although  multiplying  the  experimental  value 
for  the  C^  by  a  factor  of  two  brought  them  into  excellent  agreement, 

Fig.  14b.  It  is  not  known  why  this  is  so,  but  absolute  intensity  work  does 


-35- 


(a) 


ti 


20/isec— *|  l—  SHOCK 

FRONT 


p,  =IOtorr,  1%  CO  in  Ar 


(b) 


/i  SEC 


Fig.  9a  Experimental  data:  shock  speed  =  2.  86  mm/jjsec;  T  at  front 
=  7600°K.  Because  of  non-coincidence  of  optical  beams,  the 
traces  do  not  rise  at  the  same  part  of  the  oscillograph  sweep. 

b  Analysis  of  (a):  The  solid  lines  are  smoothed  experimental 

data.  For  CO  two  lines  are  given,  marked  A  and  B.  A  refers 
to  the  trace  height  above  the  base  line  and  B  is  obtained  making 
an  allowance  for  Bremsstrahlung  which  is  assumed  to  grow 
linearly  from  the  shock  front.  a  is  the  degreee  of  dissociation 
of  the  CO.  The  points  coded  with  open  squares,  f]  for  C?  and 
open  upright  triangle s,  A  for  CO,  always  refer  to  the  races 
given  in  Table  1.  Open  circles,  O  for  C^  and  crosses  X  for 
CO,  refer  to  a  calculation  made  with  Table  1  rates  but  ignoring 
the  boundary  layer. 
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(c) 


fj.  SEC 

Fig.  9c  A  second  analysis  of  (a):  The  steady  state  approximation  was 
used  for  the  calculation  and  boundary  layer  effects  included. 
As  before  P  and  A  refer  to  C2  and  CO  using  Table  1  rates. 

O  and  X  were  obtained  putting  k7  =  2  x  10-34>  ]<?  =  2  x  10-l0, 
k"  =  2  x  IO-33. 

d  A  third  analysis  of  (a):  The  integration  of  the  rate  equations 
allowing  for  boundary  layer  was  used.  O  and  X  were  obtained 
putting  kj  =  2  x  10” 34  ,  k?  =  4x  10" 10  ,  k"  =  2x  10‘33  , 
k"3  =  2  x  10'9  ,  k|4=2x  IQ'33  . 
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(a) 


(b) 


Fig.  10a  Experimental  data.  Shock  speed  2.  82  mm/jisec,  T  at  front  = 
7400°K.  The  B rems strahlung  is  very  evident  inthe  CO  channel 
80  jxsec  behind  the  shock  front. 

b  Analysis  of  (a)  using  rates  from  Table  1. 
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Fig.  1 


IO/i.sec-*|  |«—  SHOCK 

FRONT 


Pj  =  3  torr,  2%  in  Ar 


(b) 


a  Experimental  data.  Shock  speed  =  3.  00  mm/p.sec,  T  at  front 
=  8300°K. 

b  Analysis  of  (a). 
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Fig.  12a  Experimental  data.  Shock  speed  =  3.  94  mm/jLsec,  T  at  front 
=  7300°K 

b  Analysis  of  (a).  This  run  was  combustion  driven  and  the  derivation 
of  the  calculated  and  observed  profiles  at  some  distance  behind 
the  front  is  due  to  the  severe  attenuation  which  occurred. 
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(a) 


p(=  10  torr  ,  0.1  %  CO  in  Ar 

(b) 


Fig.  13a  Experimental  data.  Shock  speed  =  3.  08  mm/jisec,  T  at  front 
=  8800°K.  Both  traces  are  of  C£.  The  lower  is  at  lower 
amplification  than  the  upper  and  has  signals  generated  from  heat 
transfer  gauges  impressed  upon  it.  On  the  upper  trace  there  is  a 
slight  overshoot  at  the  shock  front  probably  due  to  a  trace  of 
impurity.  The  marked  rise  occurring  about  60  fj, sec  after  the 
shock  front  is  due  to  the  onset  of  ionization. 

b  Analysis  of  (a).  Three  computations  were  made.  The  open 
squares,  n.  refer  as  before  to  rates  in  Table  1.  The  open 
circles  O,  have  the  same  rates  but  kj  =  1  x  10'^  and  the 
crosses,  X,  kj  =  0.  5  x  10”34 
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p  =  15  torr  ,  0.95%  CO  f 0.05% 
CC^4  in  Ar 


(b) 
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Fig.  15a  Experimental  data.  Shock  speed  =  2.  76mm/jusec,  T  at  front 
=  7000°K.  The  channel  shows  a  strong  overshoot  which 
decays  within  4  /isec  of  the  shock  front.  This  is  due  to  the 
dissociation  of  the  CCl^  added. 

b  Analysis  of  (a).  The  open  diamonds  for  ,  and  inverted 

triangles  for  CO,  V,  were  obtained  from  the  steady  state  cal¬ 
culations  with  Table  1  vahies  except  k^  =  2  x  10"  10  and  k7 
=  2x10-33.  2  4 


-43- 


B21  10 


(a) 


20 /X  SEC 


P,sl5torr,  I  %  CO  4  1%  02  in  Ar 
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/J.  SEC 

Fig.  16a  Experimental  data.  Shock  speed  =  2.  80  mm/jusec,  T  front 
=  7200°K. 

b  Analysis  of  (a).  In  this  experiment  the  O ^  was  dissociated  within 
1  /isec  of  the  shock  front.  These  calculations  are  presented 
and  A  refer  to  Table  1  rates.  The  half  shaded  squares  □  and 
circles  O  at  the  top  have  the  same  rates  except  k  ^  =  3  x  10-10 
and  those  shaded  at  the  bottom,  □  and  O  have  kj  =  7  x  10“^®. 
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tend  to  have  scatter  of  factors  of  two  and  since  the  shape  of  the  profile  is 

good,  the  failure  to  match  the  concentration  exactly  is  not  thought  significant. 

The  possible  contribution  to  the  channel  by  CO  triplet  radiation  is 

probably  small.  It  was  estimated  from  the  photographs  of  the  spectra  that 

the  CO  radiation  if  it  occurs  in  the  visible  in  at  least  an  order  of  magnitude 

less  intense  than  the  C^  since  it  was  never  observed.  It  is  possible  that 

3  3 

there  is  a  small  component  of  triplet  band  radiation  (d  A  —  a  11  ).  The 
oscillator  strength  of  this  system  is  not  known  and  neither  are  the  mech¬ 
anisms,  if  they  exist,  for  the  production  of  the  excited  CO.  Arguments 
made  earlier  in  this  paper  indicate  that  the  most  likely  collision  partner  is 
CO  itself,  and  in  this  case  the  CO'1'  will  rise  proportional  to  [CO]  and  time 
until  about  3  /j.s  in  laboratory  coordinates.  The  computed  C^  also  rises 
proportional  to  [CO]  and  time,  once  an  initial  time  has  elapsed.  A  de¬ 
crease  in  the  [  COj  concentration  should  allow  the  effect  to  be  sorted  out, 

* 

since  the  CO  will  always  rise  in  a  time  of  order  3  ]U  sec.  Unfortunately,  the 

noise  also  increases  and  mixtures  of  less  than  l/Z%  CO  were  not  investigated. 

»•/ 

The  results  never  indicated  that  CO"  might  contribute  appreciably  to  the 
C^  channel  signal. 

The  Effect  of  Added  Oxygen 

The  analysis  shown  in  Fig.  16b  of  a  1%  CO  plus  1%  O^  plus  98%  Ar 
indicates  that  a  large  amount  of  oxygen  speeds  up  the  reaction  particularly 
near  the  shock  front.  In  this  experiment  the  O^  will  be  dissociated  in  the 
order  of  the  resolution  time  of  the  experiment.  Although  a  pathway  via  CO 
could  be  contemplated  assuming  CO  +  O  —  CO  +  O  to  be  rather  faster  than 
CO  +  CO  — ►  CO  +  CO  this  would  require  O  to  be  more  effective  than  CO 
which  already  has  a  close  to  gas  kinetic  cross  section.  This  is  improbable. 
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The  observations  are  explained  by  reaction  (3)  followed  by  (5). 

Since  the  rate  of  (5)  is  already  known,  only  the  rate  of  (3)  is  to  be  found. 

A  good  fit  is  obtained  with  kj  =  5  x  10  ^  cm^  s  \  although  the  computations 
are  not  particularly  sensitive  to  50%  changes  in  this  value,  Fig.  16b. 

Assuming  the  mechanism  of  formation  of  already  given,  a  check 
was  made  by  the  replacement  of  5%  of  the  CO  by  O^,  in  a  mixture  with  99% 
argon  at  6900°K  it  was  found  that  the  C^slope  was  a  factor  of  three  lower 
than  that  expected  from  the  relation  in  Fig.  8.  Under  these  conditions  the 
backwards  reaction  (2)  will  be  important  and  the  result  indicates  that 
k  ^ /k  3  =  2  x  10  ,  a  value  quite  consistent  with  the  rates  obtained  from 

other  expressions.  Such  experiments  further  indicate  the  appropriateness 
of  the  mechanism  given  for  formation. 

Possibility  of  CO  Being  Important 

jj( 

If  the  arguments  based  on  Ref.  15  are  invalid,  CO  could  also  be 

involved  in  the  mechanism  of  the  overall  reaction.  In  the  absence  of 

independent  rates  for  reactions  (2)  and  (4),  it  is  impossible  to  do  more  than 

* 

find  the  limits  of  rates  for  a  mechanism  based  on  CO  ,  and  for  completeness' 

sake,  computer  runs  were  also  made  assuming  that  the  incubation  was  due 
♦ 

to  CO  ,  formed  by  the  reaction 

CO  +  C  —  CO*  +  C  (13) 

and  followed  by  * 

CO  +  M  —  C  +  O  +  M  (14) 

and  it  was  re  quired  that  k~^  =  2x10^  and  k~^  =  2x10  Fig.  9d.  This 

rate  of  (13)  is  inordinately  high  and  the  true  rate  is  almost  certainly  lower. 

Further  experimental  work  will  be  necessary  before  reactions  (13)  and  (14) 
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can  be  either  eliminated  or  included  with  any  certainty.  Because  of  their 
possibility,  the  rates  found  for  (2)  and  (4)  must  also  be  considered  as  upper 
limits. 

Steady  State  Approximation 

The  steady  state  solutions  obtained  are  also  presented  in  Figs.  9c 
and  15b.  It  will  be  noted  that  they  give  quite  a  fair  representation  and  are 
not  significantly  different  from  the  full  integration  procedure  which  uses 
more  computer  time.  Although  they  have  been  in  use  for  many  years, 
steady  state  approximations  are  frequently  considered  to  be  badly  in  error. 
In  this  system  at  least,  this  is  not  true.  Its  advantage  lies  with  those  with¬ 
out  access  to  computers.  For  instance,  even  taking  wide  intervals,  the 
steady  state  method  does  not  lead  to  oscillating  solutions  and  hand  calcula¬ 
tions  are  readily  made.  The  results  are  also  presented  to  justify  the  use 
of  steady  s  ate  arguments  used  earlier.  It  is  evident  that  while  the  magni¬ 
tudes  of  reaction  rates  so  obtained  may  be  in  error,  the  approximation  gives 
the  right  form  of  the  result. 

SUMMARY  AND  CONCLUSIONS 

The  dissociation  of  carbon  monoxide  has  been  investigated  by 
observing  its  infrared  and  visible  radiation.  The  reaction  is  not  a  simple 
one  and  is  best  explained  as  a  chain  reaction  with  C^  as  the  important 
intermediate.  The  rates  which  best  fit  the  data  based  on  computed  histories 
including  boundary  layer  buildup  but  not  attenuation  are: 

Reaction  TABLE  1  Rate  Constant 

C  +  O  +  M  —  CO(X 1  2)  +  M  ,  =  2  x  10"34  cm6  s-1 

C2  +  O  —  COfX1!:)  +  C  ,  =  6  x  10“10  cm3  s-1 
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02  +  C  -*  COtX1!:)  +  O  ,  =  5  X  10"10  cm3  s"1 

C+C  +  M~C2  +  M  ,  =  3  x  10‘33  cm6  s"1 

These  rates  are  upper  bounds  which  may  need  to  be  modified  downwards 

if  alternate  reactions  prove  to  be  important.  The  uncertainties  are  estimated 

as  +  50%. 
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